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We compare the accuracy of two cluster extensions of Dynamical Mean-Field Theory in describing 
d-wave superconductors, using as a reference model a saddle-point t-J model which can be solved 
exactly in the thermodynamic limit and at the same time reasonably describes the properties of 
high-temperature superconductors. The two methods are Cellular Dynamical Mean-Field Theory, 
which is based on a real-space perspective, and Dynamical Cluster Approximation, which enforces 
a momentum-space picture by imposing periodic boundary conditions on the cluster, as opposed to 
the open boundary conditions of the first method. We consider the scaling of the methods for large 
cluster size, but we also focus on the behavior for small clusters, such as those accessible by means 
of present techniques, with particular emphasis on the geometrical structure, which is definitely a 
relevant issue in small clusters. 

PACS numbers: 71.10.Fd, 71.27.+a, 75.20.Hr, 75.10.Lp 
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Strongly Correlated electronic materials and the solu- 
tion of the models introduced to understand their be- 
havior are one of the main challenges of modern solid 
state physics. Despite the intensive activity triggered 
by the role of correlations in high-temperature supercon- 
ductors, many questions remain unanswered. The vital- 
ity of the field is testified by the development of theo- 
retical approaches designed precisely for these systems. 
Among these method, a central role is played by Dynam- 
ical Mean-Field Theory (DMFT)^, a non perturbative 
approach which generalizes the classical mean-field the- 
ory to the quantum dynamical world. The development 
of DMFT has allowed for a number of successes, start- 
ing from the first unified scenario for the longstanding 
problem of the Mott transition, and including reliable de- 
scription of the electronic properties of many correlated 
systems^. 

The idea behind DMFT is analogous to classical mean- 
field theory, namely the assumption that each lattice site 
is equivalent to any other. The difference with the static 
case is that, within DMFT, each lattice site has a com- 
pletely non-trivial dynamics. Following the above strat- 
egy, one given site can be taken as representative of the 
whole system: the lattice problem is therefore mapped 
onto a dynamical local problem, and consequently onto 
a single-impurity model. The effect of all the remaining 
lattice sites will be described by a bath, whose frequency 
dependence will be determined through a self-consistency 
condition that enforces the equivalence between the lat- 
tice and the local problems^'^. The main limitation of the 
standard (single-site) DMFT, as we have described it so 
far, is the neglect of spatial correlations. This constraint, 
indeed, introduces some limitations which are particu- 
larly relevant in low dimensionality, and in particular it 
makes it impossible to treat phases with a definite spa- 
tial ordering such as d-wave superconductivity, d-density 
waves, stripes, dimerized states. 

A few schemes have been proposed to overcome these 



limitations, re-introducing short-range correlations by 
replacing the single impurity model with a cluster- 
impurity, which contains Nc sites in a given spatial 
arrangement In this manuscript we compare two al- 
ternative schemes that represent somehow opposite per- 
spectives in their ability to describe two-dimensional cor- 
related models and d-wave superconductivity. The Dy- 
namical Cluster Approximation (DCA)^ is based on a 
momentum-space perspective and it replaces the single 
momentum-independent self-energy of DMFT with the 
set of self-energies associated to the lattice momenta of 
an iVc-site cluster. For the (cluster) impurity model, this 
approach requires periodic boundary conditions on the 
cluster. The other approach we consider, the Cellular 
Dynamical Mean-Field Theory (CDMFT), assumes in- 
stead a real space perspective, and it generalizes more 
directly the mean-field spirit of DMFT^ . In this scheme, 
a given cluster is chosen, and a "local" theory for the 
cluster degrees of freedom is obtained through the cavity 
method, replacing the effect of the rest of the lattice with 
a self-consistent effective bath. The basic approximation 
is to assume that the dynamical field experienced by the 
cluster is Gaussian. 

The properties of the two methods have been com- 
pared in several papers'" which focused mainly on the 
asymptotic behavior for large clusters, and the two meth- 
ods have been used to study many properties of the 
two-dimensional Hubbard model like, notably, d-wave 
superconductivity^. It must be underlined that, de- 
spite the simplifications introduced by the cluster meth- 
ods with respect to the full lattice problem, the cluster- 
impurity model remains a non-trivial many-body prob- 
lem, that still requires, in practice, a numerical "solver" 
in order to achieve the Green's functions. Among 
the most popular impurity solvers we remind vari- 
ous Quantum Monte Carlo methods (Hirsch-Fye deter- 
minantal method^, and the more recently introduced 
Continuous-Time Quantum Monte CarloiS), the exact- 
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diagonalization approacb^^ the numerical renormaliza- 
tion group^^. 

As a matter of fact, present "state of the art" calcula- 
tions using accurate numerical solvers are limited to fairly 
small clusters^iiH, or, if the cluster size is increased, to 
relatively small coupling and/or finite temperatur o^^'^^ , 
which may not be representative of the strong-repulsion 
regime. Such small sizes hardly allow for an accurate 
size scaling to describe the thermodynamic limit. It is 
therefore desirable to study the performance of the dif- 
ferent cluster methods as a function of the cluster size 
within some approach that allows for a solution at ar- 
bitrary values of the size. This point of view has been 
taken in Ref. 0, where a one-dimensional exactly solvable 
model has been studied using both CDMFT and DCA, 
allowing the comparison between the two methods. 

In this work we apply a similar strategy to treat a 
model which describes the essential physical ingredients 
of the cuprates, namely their two-dimensional character, 
the effects of strong correlations, and, most importantly, 
the presence of d-wave superconductivity. This model is 
the t-J model treated at a saddle point level, following, 
e.g., Ref. [l^. While this model is clearly an approxima- 
tion of the full two-dimensional t-J model (which lacks 
an exact solution), we will consider it as our "starting 
model" . In this way we will have an exactly solvable 
model, containing all the main ingredients of cuprates, 
that we can also solve using DCA and CDMFT ap- 
proaches for any finite size of the clusters. This will allow 
us on one hand to study the convergence of the methods 
in the limit of large cluster-size, but on the other hand 
it will lead to a benchmark of the methods for small and 
intermediate clusters, such as those available in present 
numerical calculation and in those that can be expected 
in a few years. Of particular interest, in this light, are the 
"geometrical" aspects of the different approaches. When 
a phase with a given spatial structure is present in a fi- 
nite cluster, we can expect different behaviors according 
to the way in which the ordered phase fits in the chosen 
cluster. This will also depend on the boundary condi- 
tions, and will mark the difference between DCA and 
CDMFT. 

The manuscript is organized as follows: In Sec. II 
we present the reference model, i.e., the saddle-point t-J 
model. In Section. Ill we introduce DCA and CDMFT, 
and their application to our model. Sec. IV presents our 
results, and Sec. V contains concluding remarks. 
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where = ^fia^af3fip is the local spin operator, 
i^i = J2afiafia IS the local clectron density and P — 
Yiii^ ~ i^i^i^il) is a projection operator which restricts 
the fermionic Hilbert space to the low-energy subspace 
of empty and singly-occupied sites; the super-exchange 
antiferromagnetic coupling J is given by /U . 

We introduce slave boson fields hi in order to keep track 
of the empty sites (holes): this representation allows, in 
fact, to replace the constraint of zero double occupancy 
with the following equality constraint. 
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where b\hi acquires the meaning of a local density of 
holes. Enforcing this constraint by means of a local La- 
grange multiplier Ai, we obtain: 
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To obtain an exactly solvable model we decouple the 
exchange interaction Si ■ Sj introducing three sets of 
Hubbard- Stratonovich fields, which allow us to treat on 
the same footing both the particle-hole and particle- 
particle channels: the reason for this approach is given by 
the requirement that the SU (2) particle-hole symmetry 
at half-filling is being preserved. 

A static mean-field approximation is then achieved by 
replacing the auxiliary fields and the Lagrange multiplier 
with their saddle-point values. Setting (S^) = and ne- 
glecting the 4-boson hole-hole interaction, which is 0{x'^) 
near half-filling (x = is the hole doping), the slave- 

boson mean-field Hamiltonian reads 



I. MODEL AND METHOD 

A. The saddle-point t-J model 

In this section we briefly review the derivation of the 
saddle-point t-J model in order to fix the notations and 
the main concepts. Even if we will not attempt to solve 
it beyond saddle point, the starting point of our analysis 
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with the particle-hole and particle-particle amplitudes 
given by 



Xij 



Ay = 



\j{fi]fjl - lilfj])- 



(5) 
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The last step of our approximation consists in decou- 
pling the kinetic term of Eq. (|4]) , and this leads to an ef- 
fective hopping amplitude t^s — t{b]jbi) for the fermionic 
degrees of freedom: assuming full boson condensation at 
T ~ 0, we can set (6]6i) = — x and consider, at 

last, teff = xt. At small doping the effects of strong cor- 
relations are thus summarized in the renormalization of 
the free-fermion hopping term, which leads to a strong 
suppression of the kinetic energy, 0{xt), and to a rel- 
ative enhancement of the super-exchange energy, 0(J), 
i.e. finite, as x ^ 0. 

In finding a self-consistent solution of the Hamiltonian 
(jll, it has been shown^^ that while at half- filling there 
are an infinite number of degenerate ground-states, con- 
nected together by the SU{2) rotations of the particle- 
hole symmetry, as soon as doping breaks the SU{2) in- 
variance the lowest-energy state is found to be^^ the d- 
wave solution Xx = Xy = X; — — A. We will 

therefore consider this kind of solution throughout our 
analysis. 

Writing the fermionic part of in Fourier space, 

we obtain 

k 

+ Ak(/^T/lkx + /-ki/kT), (7) 
Ek = -2{tcs + x) (cos + COS ky ), (8) 
Ak — 2A{cos kx — cos ky), (9) 

where the parameters x ^-nd A, representing respectively 
a renormalization of the hopping and a d-wave order pa- 
rameter, are determined by the self-consistency equations 
([SHS])- These in Fourier space read 
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A = / -— -2 (/kT/-ki)(cosfca; - cosfcy) (11) 
o J (27r) 



3 r d^fc Ak , (3Ej 



{2nf 2E^ 2 



tanh (cos k^ — cos ky 



-J 



where -Ek = y (ck ~ m) + A^ are the eigenvalues of the 
Hamiltonian ([7]). The fermionic chemical potential /i is 
instead determined by the number equation 



1-x = 1 - 



d^k (^k - m) ^^^^ /5£^k 



(12) 



(27r)' 2 
obtained by imposing the fermion density to be (l — x) 



B. Cluster Approximations 

We are now in the position to compare the exact solu- 
tion of the saddle-point model ^ with the approximate 
cluster solutions. Within cluster DMFT methods an ef- 
fective action for the cluster degrees of freedom is defined 
as 

5cff = r dr^ct,(r)g-i(r-r')^,c..(r') + 

— 3 
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where is a dynamical "Weiss" field. By computing 
the cluster Green's function Gfi^aiT) — —{Tc^cr{T)cl„), 
the cluster self-energy is obtained as 

Er(«o;„) = G-^{iuJ,,) - G-liiuj.,,). (14) 

The two methods differ in the way the new Weiss field is 
obtained through the knowledge of the cluster self-energy 

sr(*c^n). 

Within CDMFT the "local" Green's function for the 
cluster is calculated as 



1 



dk 



(15) 

where the momentum-integral extends over the reduced 
Brillouin zone associated to the iVc-site cluster, tk is 
the Fourier transform of the cluster hopping term. 
Gioc{i^nY'^ is then used to obtain a new Weiss field 



(16) 



which determines the new effective action (1131) from 
which a new G^i^{iuJn) can be obtained: the procedure 
is then iterated until convergence. We stress that this 
method does not impose lattice translational invariance. 

The spirit of DCA is instead to generalize the 
momentum-independence of the self-energy, characteris- 
tic of the single-site DMFT, to a small cluster. Thus one 
defines a coarse-grained self-energy for every reciprocal 
lattice momentum kc associated to the cluster at hand. 
The analogue of p^ . which expresses the lattice Green's 
function in terms of the cluster self-energy, is given by 



G(kc + K,iw„) = 



1 



lUJn 



UK. + K) 



(17) 

while the self-consistency relation between the Weiss field 
and the cluster self-energy now reads 



^0 ^(^c, it^n) = 



^^G(ke + K,zu;„) 



N 



K 



-HS"(kc,iw„). 

(18) 
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In these expressions, all the cluster quantities appear as 
functions of the cluster momenta k^, and the K inte- 
gration over the reduced Brillouin zone is nothing but 
a coarse-graining of the lattice Green's function around 
these momenta. 

A crucial observation is that the diagonal nature in 
momentum space of the DCA equations requires that the 
cluster part of our effective action has periodic boundary 
conditions. As we will discuss, this may represent a se- 
vere constraint, especially for small cluster sizes. 

Even if DCA is naturally defined in momentum space, 
it is useful to recover a real-space formulation also for 
this cluster approach, in order to have a unified for- 
malism which makes easier the comparison between the 
two methods. Performing a Fourier transform upon the 
cluster momenta, all the cluster quantities Q(kc) be- 
come cyclic matrices in the real-space cluster indexes, i.e., 
with the matrix elements Qij depending only on (i — j) 
modLc^ this means that translational invariance is pre- 
served within the cluster, which thus must have periodic 
boundary conditions, as we mentioned before. The real 
space formulation of the self-consistency equation p8)) is 
then given by 



'n, 



(19) 



-j^ ^ ^{iuJn + - iDCA(K) - I^^iiUn) 



K 



where [£dca(K)]^^. = ^ ^k, ^'''"'"■""^■^^(kc + K) - 
g-iK (ei-ej)^^^. ^j^^ differs from the bare i(K) used in 
CDMFT in order to satisfy the cyclicity condition, i.e., 
the translational invariance. As in CDMFT, the solution 
requires an iterative solution of a cluster-impurity model 
determined self-consistently. 

We now detail the implementation of the two ap- 
proaches for the mean-field two-dimensional t-J model. 



C. DCA for the saddle-point t-J model 

In our analysis of the t-J model, the DCA cluster 
self-energy S(kc) consists of a normal term xO^c) and 
of an anomalous term A(kc), which are respectively the 
Fourier transform of Xij and Aij within the DCA cluster; 
explicitly, they are given by 



X(kc 
A(k, 



2^DCA 



(cos kcx + COS key), (20) 
(cosfcc^ — cosfccy)- (21) 



Since our starting model is a mean-field model, these 
quantities are w-independent; however, this deficiency is 
compensated by the possibility to solve exactly both the 
lattice and the cluster problems, even for large values of 
TV,. 



The cluster parameters xS*^^ ^^"^ ^^P^ are deter- 
mined by the DCA self-consistency equations 



3 f d k 

c"^ / (/k/k)(cosfcc^ -I- cos kc 

o J (27r) 
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(cos kcx + cos key), 



8 J (27r)^ 



(/kT/-ki)(cos kc.j; - coskcy) 



^ j^J (27r)2 2EK,k^ 
X tanh (cos kc^ — cos key), 



(23) 



where eK,ke = -2tcfr [cos(fcc^ + K^) + cos{kcy + Ky)] 



X(kc) and i^K,k, = y {<^KMa - +^k,- the sec- 
ond row of Eqs. (I22H2S1) the integration over the entire 
Brillouin zone is divided into a sum over the cluster mo- 
menta and an integration over the reduced Brillouin zone. 

It is important to note that these cluster quantities, 
used in the definition of the cluster self-energy, do not 
have an immediate physical meaning and thus they can- 
not be directly compared to the corresponding lattice 
quantities of Eqs. p^ - fTTj) . The physically relevant quan- 
tities are instead given by 
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-J 
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"2" (/k/k)(cosfc^ -I- cosfcj^), (24) 



^f^u 1 4^ (/kT/-ki)(cos/c, -cosfc,), (25) 

o J (27r) 

where the expectation values must be evaluated using 
the self-consistent parameters Xd^^ and H^^p^, as in 
Eqs. ([22] [231). 

Finally, the DCA analogue of Eq. ((T^) gives the self- 
consistency equation for the chemical potential: 



1 -X 



(Pk 



iflh 



x2 \.'k 



i2n) 

^ f (fK (eK,kc - m) ^ , PE^M. 



(26) 



(27r)^ £^K,ke 



■ tanh ■ 



D. CDMFT for the saddle-point t-J model 

Considering a square cluster C with Nc = LcX sites, 
we denote hy i = i{ix,iy) = ix + {iy ^ i)Lc the cluster 
site with coordinates {ix,iy), where ix,iy — l,...,ic- 
The local cluster Hamiltonian is then obtained from the 
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lattice one by restricting all the site-index sums to the 
cluster sites: 

He ^ + ^j)flfi- - E + 

+ ^ (A^^-44+h.c.), (27) 

i,j 6 C 

where jif and t are hermitian matrices and A° a symmet- 
ric matrix. Explicitly, their expressions read 



t'ij — ^cff ^ ^ {^j,'^+V ~^ ^j.i—'f)) ' 



E 



{Xi,V^j,i+i) + Xj,fi^j,i-fi) J 



(28) 
(29) 
(30) 



Introducing the unitary matrix ?7(K) which diagonalizes 



U{K)hiK)U^iK)\ =5^,,A''(K), (35) 



we can write the cluster propagator as (the dependence 
of the matrices on K is omitted to lighten the notation) 



DciiuJn) 



IJ 



2Lt 



, (36) 



where the index / = (cr, i) denotes the (cr — l)i^ + i 
component of a Nambu spinor (cr — 1,2). 

It is now easy to get the self-consistency relations for 
the cluster parameters: 



where fj = a;, y is a lattice displacement in the x or 



y direction, 5j^i+x 



(1 X/n^»,"ic)j 



(5j,i_i(l - X]™<5i,ni<:+i) and 5j^i±y = 5j,i±L,- The self- 
consistent parameters Xi,fi and A^^jj are formally given by 
Eqs. ((SHSI), where the expectation values must be evalu- 
ated using the cluster propagator Dc, defined below. 

In order to write this propagator in a compact 
form, we first introduce the Nambu spinor 
if IV ■ ■ • ' /l^T' /U' • • ■ ! fm), which contain aU the 2Ll 
fermionic degrees of freedom within the cluster. With 
this notation, 

.,.).-(T.(.).t(0)).(^^) 

(31) 

where G„-, = ~{Tha{T)f]M) and = 
^{T fii{T)fj.^(Q)) are, respectively, the normal and 

anomalous Green's functions. We can then express Dc 
in terms of the cluster Hamiltonian parameters. 



(£)" 



iw„l — /i(K) 



(32) 



/i(K) 



i{K) + - A 
(A-)* 



(33) 

where /x is the chemical potential times the unitary ma- 
trix, t(K) is the Fourier transform of the super-lattice 
hopping matrix tR^R/, which, for |R — R'| = Lc, con- 
nects the boundary sites of neighboring clusters: 

La 

*(K)ij = tij ~ ^eff [exp(zj^ra;Lc)(5i,nL,'^j,(K-l)Le + l 

n=l 

-I- exp{-iKxLc)5j 

-I- eyiY>{-iKyLc)5i^n5j,La(La~i)+n] ■ (34) 
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/ 7 72 E i^"^. (1,0 ii,i+Lc 

(it) 
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■ ^v, {2,t+La) (2,i) 



2Li 
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[u:^(l,^)Ki2,^+l) (39) 



A,: 



g r ^ 

_3 [ <PK S 



(2,i+La 



(40) 



fmn, «<ic(ic"i) 



where f{x) = (e^ + 1)~ is the Fermi function. 

The determination of the chemical potential ^ as a 
function of the fermion density requires some more care. 
In fact, as we have seen in Eqs. ((37H40| for x and A, the 
observables are generally site-dependent in CDMFT, and 
there is no unique procedure to extract lattice observ- 
ables (translationally invariant) from cluster quantities. 
In general, for a local observable Oi, we can estimate its 
lattice counterpart with a weighted average 



Olat 



E 

iec 



(41) 



where X^iec — ^^ and we have decided to investigate 
the two extreme cases of flat average, wf**' = 1/Nc, and 
bulk value, w^""^ = 6i^b, the latter case corresponding to 
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taking the value of the observable just in the center of 
the cluster—, represented by the site b. From the local 
fermion density 



= 1- 



2Lt 



(42) 

we can therefore extract an average density rt**^' = 
1/L^^.(ni} and a bulk density n°^^^ = (rib), and we 
can adjust the chemical potential in order to satisfy ei- 
ther one of the two equations 



(flat, bulk) 



1 



(43) 



The same argument would apply in extracting the lat- 
tice self-energy parameters, to be compared with those of 
Eqs. (|10Hlip . from the corresponding cluster quantities. 
However, in considering the flat average case, we should 
note that Xij ^-i^d A^j are defined on bonds, so that for 
each direction their total number is Lc{Lc — 1) instead of 
L^; the averages will thus be given by 



flat 

A.X 



fiat 
Ay 



A fiat 
y^fiat 



' y 
y 



Xi 



i<Lc(ic-l) 



Lc{L 



1 



Lc{Lc - 1 



■ E 

i<Lc(-Lc-l) 



A, 



(44) 
(45) 

(46) 
(47) 



The bulk value estimate is instead straightforward and 
corresponds to the values of the parameters on the in- 
nermost bonds of the cluster. 



II. RESULTS 

A. Size dependence of observables 

We start our analysis of DCA and CDMFT for the 
saddle-point t-J model by considering the behavior of 
three relevant observables as a function of the linear size 
of the cluster for square lattices of Lc x Lc sites. Through- 
out this section, energy scales are expressed in units of 
J/4 and J/t — 0.4. In Fig. [T] we plot the supercon- 
ducting order parameter for different values of doping 
X = 0,0.1 and 0.2. As discussed above, for CDMFT 
one has the alternative between "bulk" and "average" 
estimates. With respect to A, we have found that the 
average estimate is much less sensitive to the size of the 
system than the bulk one, the latter being instead rather 
unreliable, compared to the exact solution, up to large 
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FIG. 1; (Color online). Normal parameter A as a function 
of the linear dimension Lc for square clusters and different 
cluster methods. Red dots are CDMFT with bulk density, 
green dots CDMFT with average density, blue squares DCA. 
The thermodynamic limit is marked by the thick green line. 
From top to bottom a; = 0, 0.1, 0.2. 



values of L^. We will thus present, in each panel, only 
average estimates of A, while considering both the aver- 
age and the bulk estimates of the electron density. The 
behavior of A as a function of the cluster size is com- 
pletely non trivial and reveals important differences be- 
tween the two methods. While DCA converges smoothly 
and faster from Lc = 5, it presents strong size-effects for 
smaller clusters: the Lc = 2 cluster overestimates A at 
large doping, while Lc = 3 and 4 produce a strongly un- 
derestimated value of A. This shows that momentum- 
space discretization is too strong to properly describe 
the spatial structure of the order parameter as long as 
the number of allowed momenta is small. On the other 
hand, while for Lc — 2 CDMFT overestimates A as well, 
the results obtained for small clusters are in general more 
reliable, with relatively small deviations from the exact 
solution and, most important, a smoother dependence on 
Lc- However, this method converges more slowly to the 
thermodynamic limit, in particular for larger dopings, 
where A is systematically underestimated. The compar- 
ison between different doping values shows indeed that 
CDMFT is quite inaccurate for x — 0.2, a doping value 
at which the hopping processes become more relevant, 
according to tcs = xt, making the system itinerant and 
consequently better described in momentum space than 
in real space. The bulk estimate of the density is typi- 
cally found to provide a better agreement with the exact 
solution. 

The same tendencies are present in the critical temper- 
ature (not shown), with a significantly enhanced overesti- 
mate of Tc compared to A in the Lc = 2 DCA cluster (we 
obtain Tc = 1.53, 1.46 and 1.30 for x = 0,0.1,0.2, while 
the exact results are 0.767, 0.601 and 0.374). The sys- 
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FIG. 2: (Color online) Normal parameter x ^ ^ function 
of the linear dimension Lc for square clusters and different 
cluster methods. Red dots are CDMFT bulk values of x, 
green dots CDMFT average values, blue squares DCA. The 
thermodynamic limit is marked by the thick green line. From 
top to bottom X = 0.1, 0.2. 



tematic underestimate of Tc for large doping in CDMFT 
is also reflected in a smaller critical doping at which su- 
perconductivity disappears {xc ~ 0.25 even for large clus- 
ters, as opposed to the thermodynamic limit Xc — 0.35). 

We finally consider the normal parameter x- Since 
this parameter coincides with A at half-filling, due to the 
particle-hole symmetry, we focus only on x = 0.1 and 0.2. 
Here, in agreement with previous studies in a one dimen- 
sional modeli, we find that in CDMFT the bulk estimate 
of the parameter is more accurate than the average one. 
Yet, CDMFT is less accurate than DCA for the cluster 
we studied, signaling that the exponential convergence of 
bulk estimates is established only for larger values of Lc- 

The different behavior between the superconducting 
and the normal parameters underlines that the accuracy 
of the different approaches, for small clusters, depends 
crucially on the quantity under consideration. In partic- 
ular, the d-wave superconducting order parameter suffers 
stronger size effects due to its peculiar structure in real- 
space, and, when the number of cluster sites becomes 
small, it is better represented by the CDMFT solution. 
On the other hand, the normal parameter, which is more 
isotropic, is better described by DCA, which favors ho- 
mogeneous states (in the sense of states without peculiar 




FIG. 3: (Color online) Doping dependence of A and x (top) 
and of the critical temperature Tc (bottom) as a function of 
doping for a 2 x 2 cluster and different cluster methods. Filled 
symbols refer to A, open symbols to x- dots are used for 
CDMFT, blue squares for DCA. The solid green line is the 
thermodynamic limit. 
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FIG. 4: (Color online) Same as Fig. [3] for 3 x 3 cluster. 



B. Doping dependence for small clusters 

In this section we focus on the smallest clusters Lc — 
2, 3, 4, 5, where we analyze in more details the doping de- 
pendence of the observables. This study is of particular 
interest because only small sizes can be handled in a full 
numerical solution of the Hubbard or t-J models using 
CDMFT and DCA approaches. For the sake of definite- 
ness, in CDMFT we use the flat average estimate for both 
A and x and the bulk estimate for the electron density. 



We start with the Lc — 2 system, namely the so-called 
2x2 plaquette. This is the smallest cluster that can 
host d-wave superconductivity, and its limited size makes 
it definitely the most popular system in CDMFT and 
DCAi^ii^. From the results shown in Fig. [3] we find that 
both methods overestimate the d-wave order parameter, 
with stronger deviations at small a; for CDMFT and for 
larger x in DCA. As far as the critical temperature is 
concerned, however, CDMFT turns out to reliably esti- 
mate the thermodynamical limit over almost the entire 
doping range, while DCA leads to a huge overestimate 
of Tc (a factor of 2 at half-filling, which increases as the 
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FIG. 5: (Color online) Same as Fig. [3] for 4 x 4 cluster. 




X 



FIG. 6: (Color online) Same as Fig. [3] for 5x5 cluster. 



doping grows). These results suggest that the geomet- 
rical constraints imposed by the 2x2 cluster have ex- 
tremely strong effects on the d-wave phase, making this 
kind of cluster hardly useful for a quantitative estimate. 
Nonetheless, the simplicity of this cluster makes it a sim- 
ple instrument to analyze the essential physics of two- 
dimensional correlated models. 

As soon as we increase the size of the cluster to Lc ~ 3, 
CDFMT experiences a substantial improvement. Both 
A and Tc are indeed reasonably close to the exact so- 
lution, except for a moderate bifurcation of Tc for large 
doping. Conversely, DCA strongly underestimates both 
quantities. It should be noted that odd values of Lc ex- 
plicitly break the particle-hole symmetry which holds at 
half-filling in the original model, since the (tt, tt) point of 
the Brillouin zone is not included among the cluster mo- 
menta: this is the reason why x ^nd A are not equal for 



a; = 0. The success of CDMFT for this small cluster has 
a twofold interest: on one hand it is a possible promis- 
ing direction for full numerical solutions of the Hubbard 
model, since the size of this cluster is reasonably small to 
allow for a reasonably accurate numerical accuracy; on 
the other hand it suggests us that, from a geometrical 
point of view, it is important to have at least two inde- 
pendent local amplitudes for the A field (the symmetry 
group of square clusters allows two independent bonds 
for Lc = 3, in CDMFT, while in the corresponding DCA 
cluster there is only one independent A(kc)). 

For the Lc — 4: cluster the two approaches give essen- 
tially analogous results, and none of them is particularly 
interesting. This underlines the fact that, for such small 
values, the precise shape of the cluster matters more than 
the total number of sites, and that the inclusion of the 
(tt, tt) point improves the DCA results, even if not by a 
large amount. 

Interestingly, at = 5 DCA becomes substantially 
more accurate than CDMFT. This confirms that DCA 
enters the asymptotic regime more rapidly than CDMFT, 
for which the finite-size effects survive to larger clusters, 
as we anticipated in the previous section. 
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FIG. 7: (Color online) Doping dependence of A and Tc in 
DCA for tilted clusters {Nc = 8, 12, 16 from top to bottom). 
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C. Tilted Clusters 

As we mentioned above, when dealing with very smaU 
clusters the DC A solution may suffer of severe size ef- 
fects, associated to the presence or absence of character- 
istic cluster momenta, of special relevance, such as (tt, tt), 
(tt, 0), (0, tt). a potential solution for this kind of sensi- 
tivity to the specific size and shape of the cluster is the 
use of specific "tilted" lattices, as shown in Ref fTsHTol . 
These clusters are compatible with the space group of 
the lattice, and at the same time are expected to dis- 
play less important size effects than the standard square 
systems. Our analysis shows that, unfortunately, the im- 
provement brought by the tilted lattices is not substan- 
tial. In Fig. El we display DCA results for = 8, 12, 16 
sites in tilted lattices, and we found that an accuracy 
comparable to that of the 5x5 square cluster is obtained 
only for Nc = 20, i.e., with almost the same number of 
sites. 

D. More specific lattices 



■» 
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FIG. 8: (Color online) Cross cluster and its embedding in the 
2-d lattice. 

One of the outcomes of our analysis so far is that, as 
long as the size of the cluster is not sufficiently large, the 
accuracy of the results is dominated by geometrical fac- 
tors. Therefore it is interesting to consider small specific 
clusters, whose size can be accessible to full numerical 
solutions, and that can minimize the geometrical frus- 
tration (or enhancement) of the rf-wave superconducting 
state. 

To this aim we studied the "cross" cluster, shown 
in Fig. [5] together with its embedding in the two- 
dimensional space, and small rectangular lattices. The 
star geometry can be considered as a good choice, since 
it can fit a d-wave "cross" of nearest-neighbor bonds. We 
find that the two approaches perform quite differently 
for this lattice. While CDMFT does not provide partic- 
ularly accurate results, DCA reproduces remarkably well 



the exact solution for both the order parameter A and 
the critical temperature Tc- Interestingly, the accuracy 
in the superconducting parameters is not accompanied 
by an equally good description of the normal self-energy 
X (results are not shown). The accuracy of DCA is as- 
sociated to the particular values of the cluster momenta 
Kc = (0,0), ±(27r/5,47r/5) and ±(47r/5, -27r/5), which 
exclude the special symmetry points (0, tt) and (tt, 0), 
but are at the same time close enough to the anti-nodal 
points to properly treat the superconducting order pa- 
rameter. On the other hand, the momenta which are 
most important to describe x not included, leading 
to a worse estimate. 

0.^1^ — ' — ' — , — ' — . — . — , — ' — ' — ' — , — . — ' — ' — , — ' — I 

■^""^^ -- T". I 

\ ^"-v b-bT CDMFT . 
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FIG. 9: (Color online) Doping dependence of A and Tc for 
the cross cluster {Nc = 5). Top panel is CDMFT, bottom 
panel DCA. 

As far as CDMFT is concerned, much closer agree- 
ment with the thermodynamic limit is reached, as shown 
in Fig. [TOl using (small) rectangular lattices. Already 
the 2x3 rectangle provides quite accurate results over 
the entire range of dopings, and the 3x4 one shows a re- 
markable agreement with the thermodynamic limit. An 
explanation for the success of such clusters could be found 
in the relatively large number of independent bonds com- 
pared to their total number, which is a consequence of 
the lower symmetry of the cluster space group. DCA on 
these rectangular lattices does not lead to particularly 
accurate results. Essentially the results (not shown) can 
be seen as a slight improvement on the corresponding 
square lattice (the largest square lattice contained in the 
rectangle), as far as Tc is concerned. 

We have thus identified at least two relatively small 
lattices which provide accurate results and that can be 
reasonably approached using a full numerical solution of 
CDMFT or DCA for the quantum Hubbard or t-J mod- 
els, namely the 5-site cross for DCA and the 6-site rect- 
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FIG. 10: (Color online) Doping dependence of A and Tc in 
CDMFT for rectangular clusters {N^ = 6 and 12 from top to 
bottom). 

angle for CDMFT. If we assume that dynamical effects 
will not spoil the geometrical effects that we have identi- 
fied, these clusters could be an ideal compromise between 
accuracy and computational effort. 

III. CONCLUSIONS 

In this paper we have investigated the accuracy of two 
cluster extensions of DMFT which are based on dia- 
metrically opposite perspectives: DCA, which enforces 
a momentum-space point of view, and CDMFT, which 
is formulated in real space. Even if the two methods 
are similar in spirit, they turn out to have different 
strengths and weaknesses. In this paper we analyzed 
the behavior of the two methods for a number of dif- 
ferent clusters, using as a reference model an exactly 
solvable one which presents d-wave superconductivity as 
well as strong-correlation effects leading to the half-filling 
Mott physics. The model under consideration is the 



two-dimensional t-J model treated within the mean-field 
slave-boson method. Even if this model is more naturally 
seen as an approximation of the real t-J model, here it is 
used to benchmark the cluster methods against an exact 
solution containing the essential physics of the cuprates. 

The choice of the model does not allow us to discuss the 
frequency dependencies of the observables. Therefore our 
strategy is to focus on the geometrical effects introduced 
by the finite-size of the clusters (and eventually by their 
shape) within the different methods. 

Analyzing the results as a function of the cluster 
size for square lattices of linear size Lc, we found that 
the smallest cluster {Lc = 2) provides rather inac- 
curate results, at least quantitatively. This suggests 
that full CDMFT and DCA studies of two-dimensional 
models could be poorly representative of the thermody- 
namic limit. The evolution increasing the cluster size is 
quite irregular, but it shows some important properties: 
CDMFT is found to adapt rather well to some precise 
shapes, while DCA shows larger geometrical effects, even 
though it converges faster to the thermodynamic limit 
(the latter is essentially reached for lattices of the order 
of 5 X 5 sites). The limitations of DCA for small clusters 
are not dramatically reduced by using "tilted" lattices, 
which include the most relevant momenta. 

On the other hand, it is found that specific small lat- 
tices can provide very accurate results. In particular rect- 
angular lattices provide rather accurate results within the 
CDMFT approach, even for the smallest case of a 2 x 3 
rectangle. As far as DCA is concerned, we find that a 5- 
sites "cross" cluster gives extremely accurate results for 
both the superconducting order parameter and the crit- 
ical temperature, even if the description of the normal 
self-energy is not equally accurate. These two clusters 
(6-site rectangle in CDMFT and 5-site cross in DCA) 
provide probably the best compromise between size of 
the cluster and computational cost, and they can be use- 
ful suggestions for future full solutions of the actual two- 
dimensional systems such as the Hubbard and the t-J 
models. 
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